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SECTION  I 


INTRODUCTION 

Surface  dust  that  has  been  scoured  and  lofted  over  a  wide  area 
surrounding  a  near-surface  nuclear  detonation  can  have  potentially 
severe,  deleterious  effects  on  communications  and  sensor  systems. 
Consequently,  the  effects  of  dust  in  nuclear  sweep-up  and  dust-pedestal 
clouds  are  incorporated  in  the  nuclear-effects  codes.  However,  many  of 
the  assumptions  concerning  both  the  mean  properties  of  the  dust  distri¬ 
bution  and  the  fluctuations  in  refractive  index  at  microwave  frequencies 
produced  by  irregularities  in  dust  density  have  not  been  verified 
experimentally,  and,  in  at  least  some  instances,  appear  to  be  erroneous 
or  inaccurate.  SRI  International,  under  Defense  Nuclear  Agency  (DNA) 
and  U.S.  Army  Ballistic  Missile  Defense  Advanced  Technology  Center 
(BMDATC)  sponsorship  (Contract  DNA  001-83-C-C203 ) ,  fielded  an  experiment 
to  measure  the  rf  effects  produced  by  the  dust  pedestal  lofted  by  the 
DIRECT  COURSE  high  explosive  detonation.  In  particular,  the  results 
from  those  measurements,  although  incomplete,  have  pointed  out  apparent 
deficiencies  in  the  application  of  conventional  turbulence  theory  to 
dust  effects. 

In  this  project,  we  put  forward  a  plan  to  complete  the  DIRECT 
COURSE  measurement  program  and  to  resolve  the  dust  effects  issues  that 
have  recently  emerged.  The  planned  work  was  designed  to  exploit  the 
MINOR  SCALE  test,  which  took  place  in  mid-June  1985.  The  planned 
microwave  experiment  was  cancelled  so  we  initiated  a  study  of  intrinsic 
mechanisms  for  toroidal  instabilities  in  the  dust  pedestal.  This  study 
was  carried  out  by  Physical  Research,  Inc.,  under  subcontract  to  SRI 


SECTION  2 


THE  MINOR  SCALE  EXPERIMENT 

The  primary  objectives  of  the  planned  experiment  were  to 

•  Make  spaced-antenna/angle-of-arrival  measurements  for  signals 
passing  along  a  number  of  paths  at  various  elevations  above  the 
surface  that  could  be  used  to  establish  the  parameters  of  dust- 
induced  signal  fluctuations  or  "turbulence"  and  the  spatial 
variation  thereof 

•  Measure  the  changes  in  signal  amplitude  and  phase  caused  by  the 
dust  loading,  then  estimate  the  magnitude  of  the  loading 

•  Complete  and  extend  the  DIRECT  COURSE  measurements. 

SRI  planned  to  make  these  measurements  with  the  rf  transmission  tech¬ 
nique  using  three  X-band  (8086-MHz)  transmitters  and  nine  receivers  at 
three  locations — thereby  providing  nine  different  paths  to  penetrate  the 
dust  cloud.  We  planned  to  modify  the  9086-MHz  part  of  the  equipment 
built  for  DNA  and  BMDATC  for  the  DIRECT  COURSE  transmission  experiment 
for  this  purpose.  Support  would  be  needed  from  FCDNA  for  testbed  con¬ 
struction  and  installation,  aircraft  tracking,  and  frequency  allocation 
approvals . 

There  were  to  be  two  phases  to  the  effort.  The  first  involved 
the  necessary  preliminary  activities  of  detailed  experimental  planning, 
design,  and  coordination,  plus  construction  of  items  (e.g.,  antenna 
supports)  to  be  delivered  to  the  test  bed  well  before  the  main  fielding 
effort.  Part  of  this  activity  included  attending  the  Project  Officers 
Meetings  scheduled  by  FCDNA  and  providing  support  to  FCDNA.  The  second 
phase  involved  the  preparation  of  the  equipment  and  the  fielding  and 
execution  of  the  experiment.  An  extensive  data-analysis  effort  was  not 
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part  of  the  program;  rather,  a  quick-look  analysis  was  proposed,  using 
some  of  the  software  developed  for  the  DIRECT  COURSE  data  analysis, 
primarily  to  assess  the  data  quality.  A  future  effort  to  reduce  the 
MINOR  SCALE  data  and  to  analyze  those  results  (incorporating  the  DIRECT 
COURSE  results  as  well)  was  to  be  proposed  following  successful 
accomplishment  of  the  experiment. 

SRI  planned  to  reconfigure  part  of  the  equipment  assembled  for  the 
DIRECT  COURSE  transmission  experiment  (under  Contract  DNA  001-83-C-0203 ) 
to  provide  for  one  airborne  and  two  ground-based  transmitters  and  nine 
receiver  channels.  Each  receiver  channel  would  accommodate  all  three 
transmitter  signals  by  means  of  frequency  multiplexing  (the  ground-based 
transmitters  would  be  slightly  offset  in  frequency  from  each  other  and 
from  the  airborne  transmitter).  Figure  1  is  a  perspective  view  of  the 
experimental  layout,  which  was  designed  to  exploit  the  precursed  dusty 
radial  on  a  noninterference  basis. 

A  formal  technical  support  plan  (TSP)  and  a  program  document  were 
submitted  to  Field  Command  DNA,  who  designated  this  experiment  as  MINOR 
SCALE  Experiment  Number  8540.  Project  personnel  attended  the  Project 
Officers  Meetings  held  during  the  weeks  of  24  September  and  26  November 
1984.  During  these  meetings,  most  of  the  details  of  the  construction 
and  logistical  support  requirements  were  finalized,  and  the  preliminary 
layout  of  the  experiment  was  established. 


Work  was  begun  on  our  task  to  reconfigure  the  equipment.  At  this 
point,  we  were  informed  that  our  experiment  would  not  be  included  in 
MINOR  SCALE,  and  that  the  remaining  effort  should  be  directed  toward 
planning  for  future  tests.  Part  of  our  subsequent  activity  was  to 
prepare  a  TSP  and  a  program  document  for  the  MISTY  PICTURE  HE  event 
scheduled  for  Spring  1986.  Both  documents  were  submitted  to  DNA.  The 
MISTY  PICTURE  event,  however,  has  now  been  delayed  until  May  1987. 

Much  of  our  work  has  been  directed  toward  designing  an  experiment 
to  address  the  issues  raised  by  the  discrepancies  between  the  results  of 
our  experiment  for  DIRECT  COURSE  and  the  theoretical  expectations  that 
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are  embodied  in  the  nuclear-effects  codes.  The  observed  turbulence 
differed  from  that  predicted  by  the  Kolmogorov  turbulence  theory,  and 
the  dust-pedestal  rise  rate  was  far  in  excess  of  that  inferred  from  the 
two-dimensional  hydrodynamic  codes  (HULL  code).  The  latter  discrepancy 
appears  to  be  caused  by  a  three-dimensional  instability  in  the  pedestal 
region.  It  is  a  significant  issue  because  it  Implies  greater  dust 
lofting  and  entrainment  than  predicted  and  is  a  mechanism  for  producing 
continued  turbulence.  As  a  result,  we  initiated  a  study  of  intrinsic 
mechanisms  for  toroidal  instabilities.  This  is  relevant  because  the 
initial  condition  before  the  instability  sets  in  appears  to  be  the  thin, 
"reverse"  torus  near  the  surface,  which  is  predicted  by  the  hydrocode 
calculation  for  1-MT  nuclear  height-of-burst  (HOB)  events.  Section  3  of 
this  report  presents  the  results  of  this  study. 


SECTION  3 


INVESTIGATION  OF  INSTABILITY  MECHANISMS 

3 . 1  BACKGROUND . 

An  in-depth  investigation  of  the  following  possible  dust-density 
irregularity  production  and  dissipation  mechanisms  was  carried  out 

•  Conditions  for  Mach  reflection 

•  Mechanism  for  vortex  formation 

•  Vortex  trajectory 

•  Vortex  instabilities 

•  Post-shock-wave  hydrodynamic  instabilities 

•  Post-shock-wave  convective  instabilities. 

This  work  reported  in  this  chapter  focuses  principally  on  vortex 
formation,  trajectories,  and  circumferential  instabilities,  which  are 
set  forth  as  the  principal  mechanism  for  formation  of  ring-wave  cell 
structure.  This  structure  was  observed  to  form  close  to  the  ground 
during  the  DIRECT  COURSE  event.  Post-shock-wave  convective  instability 
was  also  investigated  as  a  possible  mechanism  for  formation  of  the 
observed  cellular  structure. 
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3.2  VORTEX  FORMATION,  MOTION,  AND  CIRCUMFERENTIAL  STABILITY. 


3.2.1 


The  DIRECT  COURSE  Event. 


DIRECT  COURSE  was  a  600-ton  yield,  50-ra  height-of-burst  high- 
explosive  event.  The  explosive  utilized  a  mixture  of  ammonium  nitrate 
and  fuel  oil  (ANFO).  Data  recorded  during  the  event  included  photo¬ 
graphic  coverage,  dynamic  and  static  pressure,  hot-wire  measurements 
(mean  and  fluctuating  velocity  and  temperature  components),  dust-density 
levels,  and  electromagnetic-wave  propagation  characteristics. 
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In  particular,  photographic  coverage  revealed  the  formation  of  a 
"dust  pedestal"  ring  (Figure  2),  which  was  observed  to  form  close  to  the 
ground  shortly  after  detonation.  The  time  history  of  the  height  and  the 
radius  of  the  ring  for  the  first  second  after  detonation  as  derived  from 
our  data  is  shown  in  Figure  3.  Note  that  both  the  radius  and  height 
level  off  to  constant  values  of  about  200  m  and  35  m,  respectively, 
after  about  0.5  s.  In  addition,  hydrodynamic  calculations  using  the 
HULL  computer  code  reveal  a  ground-level  vortex  forming  at  about  1  s  at 
a  radial  location  of  about  180  m,  and  persisting  at  that  location 
through  times  as  late  as  15  s  (Figure  4).  Moreover,  observation  from 
above  revealed  a  circumferential  cellular  structure  in  the  dust  pedestal 
consisting  of  about  30  cells,  the  formation  of  which  is  the  main  subject 
of  the  investigation  carried  out  herein.  The  two  primary  candidates 
investigated  to  explain  the  dust  pedestal  cellular  structure  are  ring 
vortex  instabilities  and  convective  cells  (also  known  as  Benard  cells). 


3.2.2 


Mechanism  for  Vortex  Formation. 


Line  and  ring  vortices  form  as  a  result  of  shear  or  swirl  in  flows, 
which  may  be  caused  by  a  variety  of  flow  mechanisms;  e.g.,  in  the  labo¬ 
ratory  vortex  rings  are  produced  by  pulsing  a  fluid  through  a  sharp- 
edged  circular  orifice  and  the  vortex  forms  due  to  the  shearing  motion 
of  the  fluid  in  the  vicinity  of  the  orifice  edge,  which  causes  the  flow 
to  roll  up  into  a  vortex  ring.  Early-time  vortical  toroidal  motion 
within  a  rising  buoyant  thermal,  however,  develops  because  of  the 
reverse  flow  toward  the  center  of  the  thermal — as  a  result  of  inwardly- 
directed  pressure  gradients  generated  by  passage  of  the  blast  shock 
wave . 


Photograph  of  DIRECT  COURSE  event  showing  dust  pedestal  ring  15  s  after  detonation. 
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Figure  3.  Dimensions  of  the  DIRECT  COURSE  cloud  pedestal. 


One  such  mechanism  is  the  presence  of  a  "triple  point,"  which  occurs 
when  a  shock  wave  reflects  off  a  surface  and  conditions  are  such  that  a 
regular  reflection  cannot  occur,  giving  rise  to  a  so-called  Mach  reflec¬ 
tion  (Figure  5;  M  -  Mach  number,  p  -  pressure).  The  extreme  incidence 
angle,  u)^,  beyond  which  a  regular  reflection  cannot  occur,  is  shown  in 
Figure  6  as  a  function  of  shock  strength,  £  =  p^/p^,  the  reciprocal  of 
the  pressure  ratio  across  the  incident  shock  wave.  Formation  of  a 
triple  point  gives  rise  to  a  slipstream  [Figure  5(b)],  which  divides  the 
downstream  flow  above  and  below  the  triple  point.  Because  the  flows  in 
these  two  regions  traverse  different  shock  systems,  the  entropy  and 
velocity  are  different  on  the  upper  and  lower  sides  of  the  slipstream. 
The  only  quantities  that  are  constant  across  the  slipstream  are  the 
pressure  and  flow  direction.  This  velocity  difference  (the  "slip"  in 
slipstream)  occurs,  in  reality,  across  a  very  thin  region  of  high- 
velocity  gradients,  which  can  be  represented  by  a  vortex  sheet.  It  is 
‘his  region  of  high-velocity  gradients  (high  shear)  that  provides  a 
possible  formative  mechanism  for  the  ring-vortex  explanation  of  the  dust 
pedestal  phenomenon  observed  in  the  DIRECT  COURSE  event. 


90 


Figure  6.  Extreme  incidence  angle,  o;e,  beyond  which 
regular  reflections  cannot  occur  as  a  function 
of  shock  strength  £  =  p^/p2* 

To  see  at  what  radius  along  the  ground  the  blast  shock,  must  be  in 
order  for  the  triple-point  phenomenon  to  occur,  the  Incidence  angle,  id, 
is  defined  as  the  angle  that  the  tangent  to  the  expanding  spherical 
blast  shock  wave  at  its  ground  intersection  point  makes  with  the  ground 
plane  (Figure  7).  From  Figure  6  note  that,  at  early  times  (low  values  of 
shock  strength,  ?),  u  must  be  >  40  deg  in  order  for  a  Mach  reflection  to 
occur.  For  the  DIRECT  COURSE  blast  shock  centered  at  50  m  above  the 
ground,  this  angle  will  occur  at  a  radius  of  about  46  m  (Figure  8)  as 
obtained  from  HULL  code  calculations.  Physically,  this  triple  point  is 
formed  by  the  intersection  of  incident,  reflected,  and  Mach  reflection 
shocks.  Under  the  assumption  of  perfect  reflection  at  the  ground  (no 
energy  absorption),  the  shock  structure  is  equivalent  to  that  generated 
by  the  intersection  in  the  horizontal-symmetry  plane  of  the  spherical 
blast  shocks.  These  shocks  are  created  by  two  bursts  of  the  same  yield 
(detonating  simultaneously)  whose  centers  are  separated  by  a  distance 
twice  that  of  the  single  burst  above  the  ground  (Figure  9).  For  the 
DIRECT  COURSE  event,  plots  of  the  bow-shock-radlus  temporal  history 


References  are  listed  in  order  of  appearance  at  the  end  of  this  report. 
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Figure  9.  Dual  bursts  generating  shock  interaction  equivalent  to  ground  reflection 
of  single  burst. 


and  the  triple-point  trajectory  (as  calculated  by  s-Cubed  using  the  HULL 
computer  code)  are  shown  in  Figures  10  and  11.  Typical  HULL  calcula- 
tional  results  for  pressure  distribution,  from  which  bow-shock  radius 
and  triple-point  location  are  determined,  are  shown  in  Figure  12. 

To  determine  whether  the  shear  layer  generated  downstream  of  the 
triple  point  is  a  possible  causative  mechanism  for  the  ring-vortex 
explanation  of  dust-pedestal  formation,  the  velocity  difference  across 
the  triple-point  slipstream  was  considered,  as  shown  in  Figure  13,  also 
based  on  HULL  calculations.  Note  that  this  difference  is  a  maximum  at 
about  0.195  s,  which  yields  a  triple-point  radius  of  175  m,  (see  Figure 
14).  This  is  in  substantial  agreement  with  the  dust-pedestal  data 
(Figure  3),  which  show  a  radius  of  about  170  m  at  0.195  s.  Thus,  a 
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Results  of  HULL  calculations  for  DIRECT  COURSE  pressure  field  at  100,  250,  and  500  ms  showing  development 
of  triple  point. 


strong  case  is  made  for  the  triple  point  as  the  causative  mechanism  for 
ring-vortex  formation,  thereby  providing  the  basis  for  the  dust 
pedestal . 


3.2.3  Mechanism  for  Dust-Pedestal  Cellular  Structure. 

Two  possible  mechanisms  to  explain  the  observed  cellular  structure 

of  the  dust  pedestal  (Figure  15)  were  investigated,  namely,  convective 
•  » 

(Benard)  cell  formation  and  ring-vortex  instability.  Benard  cells  form 
under  conditions  wherein  a  fluid  is  heated  from  below,  causing  a  low- 
density  region  to  form  near  the  lower  surface,  giving  rise  to  convective 
instability.  HULL  code  calculations  indicate  that  the  flow  near  the 
ground  is  indeed  unstably  stratified  (density  increasing  with  distance 
from  the  ground),  which  would  support  the  Benard-cell-type  development. 
Benard  convection,  however,  falls  into  three  regimes — depending  on  the 
value  of  the  Rayleigh  number,  which  is  defined  by 


Ra 


ggATh3 

va 


(1) 


where  g  is  acceleration  caused  by  gravity,  (3  is  the  thermal  expansion 
coefficient,  AT  is  the  temperature  difference  across  the  unstably 
stratified  convective  layer,  h  is  the  layer  height,  v  the  kinematic 
viscosity,  and  a  the  thermal  diffusivity.  For  0  <  Ra  <  5600,  the 
convection  is  two-dimensional  (no  cellular  structure).  For  5600  <  Ra  < 
12,000,  the  flow  is  characterized  by  time-periodic  three-dimensional 
disturbances,  and  for  Ra  >  12,000,  more  complex  aperiodic  spatial  and 
temporal  variations  are  observed.  Indeed,  for  Ra  >  30,000,  the  convec¬ 
tion  is  sufficiently  turbulent  in  appearance  that  well-defined  thermal 
oscillations  are  rarely  discernible.  For  the  DIRECT  COURSE  event  under 
consideration,  a  typical  value  of  Rayleigh  number  based  on  a 
strat if icat ion- layer  height  of  30  ra,  a  temperature  difference  across 
the  layer  (obtained  from  HULL  calculations)  of  220  °R,  and  values  of  @, 


v,  and  a  of  0.00204/°R,  8.1  x  10  ft  /s,  and  8.3  x  10  ft  /s  respec¬ 
tively,  is  Ra  -  2  x  10^,  which  is  clearly  in  the  turbulence  aperiodic 
range,  and  thus  would  not  be  likely  to  display  the  regular  circumferen¬ 
tial  cellular  structure  observed  in  the  DIRECT  COURSE  event. 

With  regard  to  ring-vortex  instability,  the  generating  mechanism 
for  formation  of  the  ring  vortex  has  already  been  identified  in  the 
above  discussion  concerning  the  triple  point.  Periodic  waves  may 
develop  circumferentially  around  the  ring  vortex  and  have  been  experi¬ 
mentally  observed  (Figure  16).  A  theoretical  analysis  performed  by 
Widnall  and  Sullivan2  shows  the  number  of  unstable  waves  to  be  a 
function  only  of  the  ratio  of  the  radius  of  the  vortex  core,  a,  to  the 
radius  of  the  vortex  ring,  R  (Figure  17).  The  results  of  the  analysis 

are  shown  in  Figure  18,  where  n  is  the  number  of  waves  in  the  unstable 
8R 

mode  and  V  -  In  —  +  C,  where  C  is  a  constant  determined  by  the  details 
of  the  swirl  velocity  distribution.  (For  example,  if  the  velocity 
distribution  is  uniform,  C  «  -  ^).  The  abscissa  in  Figure  18,  V,  is  the 
ratio  of  the  self-induced  translational  velocity  of  the  vortex  ring,  V  . 
and  the  vortex  strength,  T,  divided  by  twice  the  perimeter,  4iR.  For 
the  30  waves  observed  in  the  DIRECT  COURSE  event,  extrapolation  of  the 
theoretical  curve  shown  in  Figure  18  results  in  V  *  5  for  n  -  30,  which, 
for  a  vortex-ring  radius  of  175  ra,  yields  a  core  radius,  a,  of  -  7.3  m. 

A  more  detailed  estimate  based  on  two-dimensional  vortex-sheet  rollup 
theory  (cf.  Section  3. 2. 3. 2  below)  yields  a  value  of  15  m.  The  core  of 
a  vortex  is  defined  as  a  region  of  constant  angular  velocity,  u>,  with 
tangential  velocity  varying  as  the  radius  according  to  v  -  ur,  where  0  < 
r  <  a  inside  the  core.  Outside  of  the  core,  the  tangential  velocity 
varies  as  1/r. 

3. 2. 3.1  Estimate  of  Vortex-Ring  Parameters  from  Triple-Point 
Calculat ions . 

A  more-detailed  estimate  of  vortex-ring  parameters,  including  core 
size,  can  be  made  utilizing  the  triple-point  data  of  Figures  11,  13,  and 
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SECTION  A-A  (Enlarged) 


14.  From  the  triple-point  calculations  of  vortex  sheet  strength  and 
vorticity  distribution,  it  is  possible  to  work  out  the  properties  of  the 
vortex  ring  formed  by  the  rollup  of  this  sheet.  Assuming  that  the  time 
scale  for  ring  rollup  is  much  larger  than  the  time  scale  for  vortex- 
sheet  formation  caused  by  shock  reflection,  the  initial  ring  radius  is 
found  by  conserving  the  vertical  impulse  (or  virtual  momentum)  in  the 
vortex  sheet.  The  effective  core  size  of  the  resulting  vortex  ring  is 
found  by  conserving  kinetic  energy.  The  total  circulation  is  found  by 
integrating  the  strength  of  the  vortex  sheet.  The  circulation  and  the 
effective  core  size,  together  with  the  ring  radius,  determine  the  self- 


induced  motion  of  the  ring  toward  the  ground.  The  circulation  and 


distance  above  the  ground  determine  the  velocity  along  the  ground  that 
the  ring  would  have  without  the  inflow  accompanying  the  explosion. 


Estimates  can  then  be  made  of  the  spatial  wavelengths  for  the  vortex¬ 


ring  instability  and  the  time  scale  over  which  it  would  occur. 


To  obtain  an  estimate  of  vortex  properties,  the  numerical  results 


for  the  strength  of  the  vortex  sheet  vs.  time,  as  well  as  the  location 


of  the  triple  point  with  time,  were  fit  with  simple  mathematical  func¬ 
tions  so  that  the  spatial  distribution  of  vorticity  y(r,y),  where  y  is 


the  vertical  direction,  could  be  obtained.  From  curve  fitting  to 


Figures  11,  13,  and  14,  it  was  determined  that  the  formulae 


h(r )  -  20(r/100)‘ 


y(t)  -  300(t-0.05)(0.525-t)/t 


r(t)  -  325(t-0.045)' 


r- 1^.' 


gave  a  reasonable  fit  to  the  numerical  results  consistent  with  the  goals 
to  provide  a  simple  model  of  vortex  behavior.  With  these  formulae, 
y(r,y)  was  available,  as  well  as  the  shape  of  the  sheet,  h(r). 

The  geometric  configuration  of  the  initial  vortex  sheets  is  shown 
in  Figure  19.  The  sheet  begins  at  a  radius  of  55  m  and  extends  to  a 
radius  of  255  m  in  meridional  (h,r)  planes.  The  height  above  the  ground 
varies  nonlinearly  with  r  in  Eq.  (2);  thus,  the  vortex  sheet  is  not 
flat.  To  simplify  the  calculations,  this  sheet  was  assumed  to  be  flat 
and  to  make  an  angle  of  32  deg  with  the  ground.  The  velocity  difference 
across  the  sheet,  Au  from  Figure  13,  was  taken  as  the  strength  per  unit 
length  of  the  vortex  sheet.  The  calculations  were  done  in  the  r'  coor¬ 
dinate  system  along  the  "flat"  vortex  sheet  (see  Figure  19). 


3. 2. 3. 2  Equivalent  Core  Size. 

Calculations  of  the  equivalent  core  size  that  result  from  the 
rollup  of  the  vortex  sheet  are  made  by  equating  the  kinetic  energy  in 
the  sheet  configuration  to  that  of  an  equivalent  vortex  core  of  constant 
vorticity.  This  would  give  (for  example)  the  correct  core  size  to  use 


if  the  self-induced  velocity,  v^,  of  the  ring  were  to  be  calculated  by 
Kelvin's  formula 


v.  -  r [ ln(8R/a)  -  1/4]/4tR  ,  (5) 

where  T  is  the  vortex  strength,  R  the  ring  radius,  and  a  the  effective 
core  size.  This  core  size  can  also  be  used  to  give  an  estimate  of  the 
wavelength  of  the  vortex-ring  instability,  although  the  calculation  is 
not  precise  because  the  details  of  the  vortex  ring  instability  depend 
upon  the  actual  distribution  of  the  vorticity  in  the  core3;  this  vortex 
sheet  will  not  roll  up  to  form  a  vortex  of  constant  vorticity.  These 
issues  are  discussed  in  Reference  4.  The  work  of  Widnall,  et  al5  shows 
that  effective  core  size  is  determined  from  the  kinetic-energy  content 
of  the  flow  in  relation  to  an  equivalent  vortex  of  constant  vorticity. 
The  configuration  for  the  calculation  of  kinetic  energy  is  shown  in 


Figure  20.  In  general,  the  kinetic  energy  per  unit  mass  of  a  flow 
containing  a  vortex  sheet  is  given  by 

KE  =  -J  A$vndS  ,  (6) 

where  A#  is  the  Jump  in  velocity  potential  across  a  sheet,  and  v  is  the 

n 

local  velocity  normal  to  the  sheet,  if  the  total  vorte*  strength  in  the 
sheet,  r,  is  not  zero,  the  contour  s  must  be  taken  to  infinity  because 

A$  does  not  go  to  zero  for  such  a  flow  at  the  edge  of  the  sheet. 

For  a  vortex  core  of  constant  vorticity,  the  kinetic  energy  of  the 
vortex  core  must  be  included.  Therefore,  for  the  kinetic  energy  per 
unit  mass  of  a  vortex  sheet  we  have 

KE  =  -/A*vndr'  .  (7) 

while  for  the  kinetic  energy  per  unit  mass  of  a  vortex  core  of  constant 
vorticity  we  have 

KE  .  -I>vndr'  *  v#2d*  .  (8) 

where  v0  is  the  tangential  velocity  and  A  is  the  cross-sectional  area  of 
the  vortex  core.  Because  the  potential  Jump  across  the  sheet  for  an 
isolated  vortex  is  constant  and  equal  to  -r,  the  integral  may  be 
evaluated  analytically.  For  a  core  velocity  of 

v(r ' )  =  rr'/2ia2  ,  (9) 

the  contribution  to  the  kinetic  energy  per  unit  mass  from  the  core  is 


y  f 


so  that  the  total  kinetic  energy  per  unit  mass  of  an  isolated  vortex 
core  of  constant  vorticity  is 

KE  =  1i'“  r2(ln  r’-ln  a  )/4i  +  r2/16i  .  (11) 

r 

Note  that  this  expression  does  not  converge  because  In  r’  continues  to 
increase  as  r'-»«.  However,  the  kinetic-energy  integral  for  the  vortex 
sheet  also  contains  the  same  In  r'  term;  when  the  kinetic  energy  per 
unit  mass  of  the  two  equivalent  flows  are  equated  these  terms  will 
cancel . 

Equating  the  kinetic  energy  per  unit  mass  from  the  vortex  sheet  to 
the  kinetic  energy  per  unit  mass  of  an  equivalent  vortex  core  of 
constant  vorticity  gives,  after  some  manipulation,  an  expression  for  the 
core  size,  a,  of  the  form 

In  a  =  -[4»KE/r2-ln  r']  +  1/4  .  (12) 


We  can  now  let  r'*«  and  obtain  a  finite  value  for  In  a.  Actually,  the 
calculations  were  done  numerically,  and  were  terminated  when  the  term  in 
brackets  reached  a  constant  value.  When  this  was  done  for  the  vortex 


sheet  in  question,  values  for  the  vortex-core  size  of  a  =  42  m  and  a 


2 

circulation,  r  (where  r  =  J  .  vdr*),  of  about  9600  m  /s  were  obtained. 

J  QnPPr 


sheet1 

This  assumes  that  the  vortex  sheet  was  created  near  the  ground  and  was 
allowed  to  roll  up  in  its  entirety  without  any  other  influences.  It  is 
clear,  however,  that  the  vortex  sheet  was  subject  to  strong  inflow  from 
the  flow  rushing  toward  the  axis.  Therefore,  it  is  unlikely  that  the 
entire  vortex  sheet  was  able  to  roll  up  as  a  unit.  If  we  assume  that 
only  the  25-m  extent  of  vortex  sheet  on  either  side  of  the  vorticity 


maximum  at  r  =  175  m  rolls  up  as  a  single  vortex,  values  of  a  =  15  m  and 

2 


vortex  strength  (circulation)  of  4330  m  /s  are  obtained.  Based  on  the 
DIRECT  COURSE  data,  there  is  no  way  to  accurately  determine  an 
experimental  value  of  a.  However,  a  value  of  15  m  compared  with  the 


overall  effective  vortex  cross-sectional  radius,  Rq  (Figure  17),  of  the 
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dust  pedestal  of  about  25  m  as  determined  from  experimental  data  (e.g., 
Figure  15)  appears  reasonable. 


When  these  parameters  are  used  in  the  self-induced  velocity 
formula,  a  vertical  speed  of  8.23  ra/s  is  predicted,  which  is  certainly 
compatible  with  the  result  shown  in  Figure  4  (HULL  calculations)  where 
the  vortex  is  shown  to  have  descended  from  an  altitude  of  50  m  at 
t  *  1.5  s  to  an  altitude  of  20  m  at  t  -  5  s. 

3. 2. 3. 3  Stability  Issues. 

As  discussed  in  Reference  3,  the  vortex-core  instability  depends 
upon  the  details  of  the  distribution  of  vorticity.  For  constant 
vorticity,  the  wavenumber,  k  (s  2t/X,  where  X  is  the  wavelength)  at 
which  the  instability  will  occur  is  given  by  ka  -  2.4;  for  distributions 
that  have  a  vorticity  maximum  at  the  center  of  the  vortex,  values  of  ka 
of  1.7  to  3.8  have  been  obtained.4  The  experimental  observations 
discussed  here  give  an  instability  having  30  waves  at  a  ring  radius  of 
about  180  m,  yielding  a  value  of  k  -  30/180  ra  -  0.167/m.  For  a  vortex 
core  size  of  15  m,  this  gives  a  value  of  ka  of  2.5 — very  much  in  the 
expected  range. 

3. 2. 3. 4  Amplification  Rate. 

An  isolated  vortex  ring  has  a  circumferential  wave  instability 

2 

amplification  rate  proportional  to  T/(4i  R  ),  which  is  proportional  to 

the  local  strain  caused  by  the  ring  flow  field.4  However,  the  principal 

source  of  instability  in  the  experimental  situation  comes  from  the 

vortex  proximity  to  the  ground;  in  this  case,  the  vortex  feels  the 

strain  field  from  its  image,  more  like  the  short-wave  vortex-pair 

instability.7  Thus,  the  amplification  rate  will  be  proportional  to 
2 

r/[2*(4h  )]  because  the  vortex  is  located  a  distance  2h  from  its  image. 
Estimating  the  height  of  the  vortex  at  20  ra  above  the  ground  gives  an 
amplification  rate,  a,  of  0.43  s  Experimental  observations2  suggest 
that  the  instability  is  complete  at  about  at  -  2 ,  or  at  a  time  of  about 


m 


•a 


m 


01 

m 


4  s  after  vortex  formation  is  complete.  Short-wave  vortex  instabili¬ 


ties,  whether  caused  by  the  strain  field  of  the  ring  or  by  the  proximity 
to  the  ground,  occur  at  the  same  wavenumber,  namely  that  of  the  critical 
second  bending  mode  as  discussed  in  Reference  3.  For  the  DIRECT  COURSE 
event,  HULL  calculations  show  the  vortex  to  be  completely  formed  at 
about  1  s  (Figure  4)  and  the  overhead  photograph  of  Figure  15  shows  the 
circumferential  cellular  structure  to  be  completely  formed  by  about 
4.5  s,  thereby  substantially  agreeing  with  the  laboratory  results  of 


Reference  2. 


3. 2. 3. 5  Impulse — Initial  Ring  Radius. 


The  vortex  rollup  method  can  also  be  applied  to  determine  the 


initial  radius  of  the  ring.  This  requires  that  we  ignore  both  the 
influence  of  the  ground  and  the  strong  inflow  upon  the  rollup.  A  vortex 
sheet  rolling  up  conserves  its  impulse  or  virtual  momentum.  The  impulse 


of  a  circular  vortex  sheet  is 


I  -  2i  / A4>  r  dr' 


(see  Figure  19)  for  definitions.  For  this  ring,  we  require  that  the 
initial  impulse  equal  the  final  impulse  of  the  rolled-up  ring,  which  is 


simply 


i  *  ri  r 


where  R  is  the  final  radius  of  the  ring.  Integrating  the  impulse  of  the 


vortex  sheet,  equating  it  to  the  impulse  of  the  ring  in  its  final  con¬ 


figuration,  and  solving  for  R  gives  R  -  161  ra,  which  is  in  reasonable 


agreement  with  the  HULL  calculations  that  show  a  vortex  at  a  radius 


180  m. 


3. 2. 3. 6  Interpretation  of  HULL  Calculations  in  Terms  of  a 


Vortex  Model . 

The  vortex  model  presented  here  is  somewhat  idealized  in  comparison 
with  the  complexity  of  the  flow  predicted  by  the  HULL  calculations  and 
shown  in  Figure  4;  it  cannot  be  expected  to  predict  all  of  the  observed 
phenomena.  Nevertheless,  it  seems  to  have  the  predictive  power  to 
explain  features  of  both  HULL  calculations  and  experiments.  Additional 
observations  can  be  made  about  the  numerical  results  based  on  vortex 
modeling.  First,  the  time  from  t  -  1  to  5  s  seems  to  be  a  time  of 
vortex  rollup.  However,  the  vortex  sheet  is  subject  to  a  strong  inflow 
toward  the  axis,  which  is  not  in  the  simple  model  and  may  affect  the 
portion  of  the  vortex  sheet  that  forms  the  concentrated  core.  The 
period  from  t  -  5  to  t  -  10  s  seems  to  be  a  period  of  vortex  motion  and 
instability.  If  the  vortex  were  isolated  in  the  presence  of  the  ground 
plane,  it  would  move  out  along  the  wall  at  a  velocity  given  by 


ground 


-  T/(4»h) 


or  at  about  17  m/s.  That  the  vortex  remains  stationary  at  r  -  180  m, 
from  t  -  5  to  t  -  15  s,  means  that  the  inflow  must  have  a  value  close  to 
this,  because  the  vortex  is  convected  by  the  combination  of  self-induced 
and  inflow  velocity. 


3. 2. 3. 7  Direct  Calculations  of  Vortex  Sheet  Rollur 


A  direct  simulation  of  the  rollup  of  the  vortex  sheet  created  by 
the  triple-point  shock  reflection  can  be  used  to  verify  the  model  of  the 
previous  section  and  to  compare  with  the  HULL  calculations.  To  obtain  a 
simple  model,  a  local  two-dimensional  vortex  sheet  was  used.  This  model 
does  not  include  the  self-induced  motion  caused  by  curvature,  but  it 
should  be  accurate  for  a  preliminary  investigation.  The  effect  of 
proximity  to  the  ground  is  included  through  a  system  of  images.  Since 
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the  velocity  field  associated  with  the  curvature  of  the  vortex  ring  is 
not  included,  and  this  is  negative  at  the  vortex  sheet,  the  height  of 
the  sheet  above  the  ground  will  be  overestimated. 


The  vortex  sheet  is  represented  by  a  number  of  discrete  vortex 
elements;  in  this  case,  100  elements  were  used.  The  velocity  of  each 
element  is  given  by  the  influence  of  the  other  vortices  including  the 
images.  The  effect  of  inflow  toward  the  axis  is  included  by  adding  a 
negative  horizontal  velocity  of  15  m/s  to  each  vortex  element.  (This 
value  was  chosen  somewhat  arbitrarily;  therefore,  the  value  of  r  at 
which  the  vortex  sheet  is  found  is  not  accurate  in  time.)  The  vortex 
points  are  moved  in  time  with  this  combined  velocity.  Proceeding  in 
time  produces  a  simulation  of  the  behavior  of  the  vortex  sheet. 


K* 


If  the  vortex  sheet  is  represented  as  a  set  of  point  vortices  of 
strength,  y^,  located  at  the  point,  (r^,  y^),  the  velocity  components  in 


the  r,  y  directions  of  each  point  are  given  by 


dr . 


— i  -  -  L  rY 

dt  2t  f-Tj 


y.  -  y. 

_ _J 


y.  +  y  . 


(yi"  yj)  +(IVrj 


+  dr)‘ 


(yt  +  y^2  ♦  (|r.-  r.|  +  dr)2 


(16) 


dy. 


r  *  57  pj  (ri-  y 


- - - - - 2 - 5 - - - o)  (17) 

(yL-  yj)  +  <1^-  rj  I  +  dr)  (y.+  y^)  +  (|r.-  r .  |  +  dr) 


where  dr  is  a  convergence  factor  that  can  be  interpreted  as  an  effective 
core  diameter  for  the  vortex  element;  this  is  a  standard  technique  in 
vortex  methods  to  avoid  chaotic  behavior. 
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Equations  (16)  and  (17)  were  integrated  numerically  from  t  -  0  to 
t  -  8  s.  The  results  are  shown  in  Figures  21  to  26.  Figures  21,  22, 
and  23  show  the  position  of  the  vortex  sheet  as  a  function  of  time;  t  - 
0  corresponds  to  the  initial  vortex  sheet  position;  the  position  of  the 
sheet  at  subsequent  times  is  also  shown.  The  formation  of  a  vortex  is 
noticeable  at  t  -  3  s.  One  interesting  feature  of  the  rollup  of  this 
vortex  sheet  is  the  strong  interaction  with  the  ground.  The  vorticity 
that  rolls  up  to  form  the  vortex  comes  from  the  left  edge  of  the  orig¬ 
inal  vortex  sheet  after  its  strong  interaction  with  the  ground.  The 
central  portion  of  the  sheet  does  not  roll  up  into  the  vortex  within 
8  s.  The  calculations  break  down  after  some  time  in  the  region  of  the 
vortex  core;  more  sophisticated  calculations  would  need  to  be  done  to 
overcome  these  difficulties. 

The  velocity  field  from  this  vortex  sheet  is  shown  in  Figures  24, 
25,  and  26  at  t  -  0,  3,  and  8  s.  The  shape  of  the  vortex  sheet  is  also 
shown  in  these  figures.  Recall  from  Figures  13  and  14  that  the  maximum 
of  the  vorticity  occurs  at  about  r  -  175  m.  This  results  in  a  strong 
vortex  flow  even  when  the  sheet  is  not  formed.  At  later  times,  the 
concentrated  vortex  rollup  becomes  a  center  of  vorticity. 

3. 2. 3.8  Conclusions . 

These  vortex  sheet  calculations  demonstrate  the  process  of  rollup 
that  forms  the  concentrated  vortex  from  the  initial  vortex  sheet  caused 
by  the  triple-point  shock  reflection.  Although  not  all  of  the  observed 
features  of  the  flow  field  are  reproduced  in  these  simple  calculations 
(because  several  important  flow-field  features  were  not  included  in  the 
model),  one  important  feature  is  clear:  not  all  of  the  vortex  sheet  is 
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Figure  23.  Calculated  vortex  position  as  function  of  time,  t  =  0  to  8  s. 
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Figure  24.  Vortex  sheet  velocity  field,  t  =  0  s. 
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Figure  25.  Vortex  sheet  velocity  field,  t  = 
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Figure  26.  Vortex  sheet  velocity  field,  t  5  8  s, 
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rolled  up  to  fora  the  concentrated  core  near  the  ground.  Therefore,  its 
circulation  will  not  be  the  circulation  of  the  initial  sheet,  and  its 
core  size  will  be  only  a  fraction  of  the  core  size  that  would  result 
from  rollup  of  the  entire  sheet.  More  detailed  calculations  could  be 
done  to  overcome  the  limitations  of  two-dimensional  flow  and  to 
incorporate  numerical  techniques  that  would  allow  the  rollup  to  be 
followed  for  longer  times. 

Thus,  when  one  considers  both  the  Benard  convective-cell  phenomenon 
and  vortex-sheet  rollup  with  vortex-ring  instability,  the  latter 
phenomenon  appear  to  be  the  most  likely  explanation  for  the  observed 
DIRECT  COURSE  dust  pedestal  and  its  ring-wave  cellular  structure. 
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